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O \ Abstract 

We investigate the evolution of the electroweak phase transition, 
^ ■ using a one-Higgs effective potential that can be regarded as an ap- 

\ proximation for the Minimal Super symmetric Standard Model. The 

phase transition occurs in a small interval around a temperature Tt 
below the critical one. We calculate this temperature as a function of 
the parameters of the potential and of a damping coefficient related 
$_i ' to the viscosity of the plasma. The parameters that are relevant for 



baryogenesis, such as the velocity and thickness of the walls of bubbles 
and the value of the Higgs field inside them, change significantly in 
the range of temperatures where the first-order phase transition can 
occur. However, we find that in the likely interval for Tt there is no 
significant variation of these parameters. Furthermore, the tempera- 
ture Tt is in general not far below the temperature at which bubbles 
begin to nucleate. 

1 Introduction 

The electroweak phase transition is an appealing scenario for the generation 
of the baryon asymmetry of the universe (BAU), since it contains the three 
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necessary ingredients known as Sakharov's conditions, namely, baryon num- 
ber violation, C and CP violation, and a departure from thermal equilibrium. 
However, to obtain a quantitatively successful electroweak baryogenesis, one 
has to consider some extension of the Standard Model (SM) in which there is 
enough CP violation as well as a sufficiently strong first-order phase transi- 
tion (see |IJ for reviews). The order parameter used to measure the strength 
of the elecroweak phase transition is the mean value of the Higgs field in the 
broken symmetry phase at the critical temperature, v (T c ) = {4>)t c - 

In the standard mechanism for electroweak baryogenesis, which assumes 
a first-order phase transition, the non-equilibrium condition acts in two dif- 
ferent ways. The first of them is through the expansion of bubbles of the 
stable phase, which combined with CP violation inside the walls of bubbles 
produce non-equilibrium particle densities, which then give rise to the baryon 
asymmetry. The second effect of the departure from equilibrium is connected 
to the baryon number violation condition: baryon number violating processes 
must be turned off before the system reaches thermal equilibrium in order to 
avoid the washout of the generated BAU. 

In the expansion of bubbles an asymmetry between left handed quarks 
and their antiparticles (and an opposite right handed asymmetry) is built up 
in front of the moving walls. This asymmetry biases the anomalous, baryon 
number violating sphaleron processes present in the high temperature phase. 
Thus, bubble walls must have a nonvanishing velocity to generate a net 
baryon asymmetry. On the other hand, the left-handed asymmetry injected 
in front of the wall needs some time to diffuse and bias the sphaleron processes 
before the reflected particles are re-caught by the wall. If the wall moves too 
fast, there won't be enough time for sphalerons to produce baryons. As 
a consequence, the generated baryon asymmetry has a peak at some small 
velocity which depends on the time scales associated to particle diffusion and 
baryon number violation p], || || . In particular, it was shown in Ref . |4[] that 
for the MSSM the peak is at v w = 0.01 - 0.03. 

The chiral quark asymmetry generated in front of the bubble wall depends 
also on the wall thickness. If the wall width is small compared to the inverse 
temperature, the CP- violating reflection of particles by the bubble wall can 
be treated quantum mechanically f|. Otherwise, the effect of CP violation 
acts as a classical force on the quarks as they pass through the wall . This 
is the case in the electroweak phase transition, since the wall width has been 
estimated to be / w ~ lOT ^ [^, ^|, Q . The final baryon asymmetry will thus 
also depend on l w , resulting a larger asymmetry for thinner walls. 

Regarding the avoiding of washout, sphaleron processes must be sup- 
pressed in the broken phase so that baryons that are created in front of the 
moving wall are not erased after entering the bubble. This requirement im- 
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poses a constraint on the sphaleron energy which can be expressed as the 
well known condition for the Higgs mean value [10 



^>1. (1) 

This condition gives a severe constraint for the theory of the electroweak 
interactions. Indeed, the minimal SM is unable to explain the observed BAU, 
since in this model the phase transition is not of the first order, but only 
a smooth crossover [fLl]| . Furthermore, electroweak baryogenesis constrains 



the Minimal Supersymmetric Standard Model (MSSM) to a small region of 
parameter space, where the Higgs mass is less than about 105GeV and the 
stop mass is less than the top one [12|, [Tj. 



From the above it becomes apparent that phase transition dynamics plays 
a relevant role in electroweak baryogenesis. Being the baryon production very 
sensitive to the bubble wall width and velocity, and to the Higgs mean value, 
it is important to determine as accurately as possible the values of these 
parameters at the actual temperature of the phase transition. Indeed, all 
of these parameters vary significantly in the temperature interval in which 
the first-order phase transition can take place, that is to say, between the 
critical temperature T c at which the two phases have the same free energy 
and the temperature To at which the barrier between the two minima of free 
energy disappears. The exact temperature T t of the transition depends on 
the evolution of the bubbles after they are nucleated, which in turn depends 
on the viscosity of the plasma. 

In this paper we will analyze in detail the dependence on temperature 
of the parameters l w , v w , and v (T) /T, and we will determine the relative 
position of the temperature T f in the interval Tq—T c . We will do this compu- 
tation for different values of the parameters of our model, placing emphasis 
on those that adequate to the context of the MSSM. We will also discuss sev- 
eral aspects of the phase transition dynamics, such as transitory states in the 
evolution of the Higgs field, or the possibility that the transition occurs near 
the temperature T , where the mechanism of baryogenesis would be different 
from the usual one. 

The plan is the following. In the next section we introduce our model 
for the first-order elecroweak phase transition, and calculate the shape of the 
nucleated bubbles as a function of temperature. In section [3] we study the 
subsequent evolution of the bubbles in the hot plasma. In section f| we deter- 
mine a temperature interval outside which the phase transition cannot occur, 
independently of the bubble evolution. Finally, in section [5] we compute the 
temperature and duration of the transition and discuss the implications for 
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baryogenesis. In the appendix we show some details of the numerical calcu- 
lation. 



2 The phase transition 

The electroweak phase transition takes place when the expectation value of 
the Higgs field passes from its high temperature value ((f)) = to its nonzero 
value in the low temperature broken phase. We will use a simple model for 
the phase transition, which is nevertheless representative of the electroweak 
theory at high temperature. 



2.1 The effective potential 



We will assume that the free energy has the well known general form & O, WA 



V (0, T) = D (T 2 - T 2 ) 2 - ET(f) 3 + A0 4 , (2) 

which will be suitable for our analysis, since it contains the essential features 
of the first order phase transition. Moreover, it can be a very good approxi- 
mation to the actual effective potential. This is for instance the form of the 



perturbative high temperature effective potential in the SM j8], p~3| . Even in 
the MSSM, which has two Higgs doublets, the free energy takes the SM-like 
form (El) in the limit in which the pseudoscalar particle of the Higgs sector is 



heavy (m^ ^> T c ) [14| 



All the parameters in Eq. (Q) depend on the particle content of the theory. 
Parameter D contains contributions from all the particles that acquire their 
masses through the Higgs mechanism. These contributions are of the form 
m 2 /v 2 , where is the zero-temperature mass of particle i and v is the zero- 
temperature vacuum expectation value of the Higgs field. Parameter E has 
only boson contributions, of the form m 3 /v 3 . In the SM we have D ~ 10 -1 , 
E ~ 10~ 2 , while in the MSSM, due to the larger particle zoo, D and E can 
be more than an order of magnitude larger than in the SM. Parameter A is 
in general temperature dependent, but it is almost constant in the range of 
temperatures in which the phase transition can take place. This parameter 
is very sensitive to the Higgs mass and for our discussion we will assume it 
to be given parametrically by m# ~ \v 2 . 

The cubic term in V ((f), T) is responsible for the first-order feature of 
the phase transition, by causing the coexistence of two minima separated 
by a barrier. Hence, the strength of the transition depends on the value of 
parameter E. At high temperature the global minimum of the potential is 
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at <p = 0. At the critical temperature 

T c = (3) 
'■t _ 

1 XD 

the two minima become degenerate, and below this temperature the stable 
minimum of V is at 




At temperature To the barrier between minima disappears, and = be- 
comes a maximum of the potential. The exact value of To depends on the 
parameters of the theory, but we can assume it to be roughly ~ lOOGeV^ 
since the dynamics of the phase transition is not very sensitive to the value 
of this constant. The number E 2 /XD is in general small, and the difference 
between T c and T is AT < 1CT 2 T C . However, as we shall see, things change 
rapidly as the temperature falls from T c to T . Hence it proves useful to 
define a dimensionless variable 

= T c 2 - T 2 T c -T 
£ ~T?-T^T C -T ' [b) 

which goes from to 1 as T runs between T c and To. 

At the critical temperature, Eq. @ gives v (T c ) = 2ET C / A. Inserting this 
value into the condition for avoiding the washout of the baryon asymmetry 
[Eq. (fj)], gives the constraint A < 2E. In the SM this results in the bound 
on the Higgs mass ttih < 40GeV, well below the experimental bound m e ^ p > 
95GeV. In the MSSM, the one loop stop contribution increases the value 
of parameter E, but this contribution is limited by the danger of inducing 
charge and color breaking minima. As a consequence, the bound on can 
be shifted only to ~ 80GeV, provided that the light stop mass is in the range 



lOOGeV <rrii<m t [0. It turns out that two loop corrections produce an 
enhancement of the strength of the phase transition, leading to an upper 
bound on the Higgs mass of ~ 105GeV |T2|. |7J . 

Although Eq. (||) is written in terms of the critical temperature T c , this 
condition is somewhat inexact, since at T c the nucleation of bubbles has not 
yet begun. On the other hand, the phase transition must be completed when 
the temperature of the Universe is still above T , because bubbles of the true 
phase are copiously produced as the barrier between minima gets small and 
disappears. Hence, the transition takes place at an intermediate temperature 
T < T t < T c with 

v(Tc) viJd <7o) ( , 

T, T t Ti 1 ' 
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The temperature T t is usually estimated by the condition 5*3 (T t ) /T t ~ 
130 — 140, where S3 is the fluctuation in free energy which is necessary for 
bubble formation H, [TBI III], 1X51, IHI. However, as was stressed in Ref. fllS 



a more careful determination of the temperature T t is important, because 
the rate of anomalous baryon number violation is an exponentially sensitive 
function of v (T) /T. In the effective potential (|2]), v (T) /T varies from 2E/X 
at temperature T c to 3E/X at To. This variation by a factor of 3/2 gives an 
uncertainty of a 50% in the bound on the Higgs mass coming from Eq. (|1]). 
Elimination of this uncertainty is thus essential given the current experimen- 
tal bounds on m#. So the correct condition for avoiding the washout of the 
baryon asymmetry would be, according to Eq. (|), 

\<p(T)E, (7) 

where p (T) is a number between 2 and 3, to be evaluated at the temperature 
of the phase transition, T t . 



2.2 Bubble nucleation 



At a temperature below T c , bubbles of the low temperature phase begin to 
nucleate. The thermal tunneling probability (per unit volume and time) for 
bubble nucleation is 



17 



r ~ A (T) e~ Ss/T . (8) 

The prefactor A (T) is roughly of order T 4 (so we will set A (T) = T 4 ), 
and S3 (T) is the three-dimensional instanton action, which coincides with 
the free energy of the nucleated bubble. At high temperature the bounce 
solution is O (3) symmetric, and the corresponding action takes the simple 
form 

r i 



r dr 



+ V 



r 



T) 



2\dr 

where r 2 = x 2 , and the equation for the configuration of the bubble is 



(9) 



d 2 (p 2 d<p 
dr 2 r dr 



V{4>) 



(10) 



with the boundary conditions (r = 00) 



and f 

dr 



0) = 0. 



For temperatures close to T c , the width of the bubble wall at the moment 
of formation is much smaller than its radius, and S3 can be expressed as a 
function of the bubble radius r , the energy difference AV between the two 
minima of the potential, and the bubble wall surface energy o [O. The 
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radius of the critical bubble to be nucleated can thus be obtained by finding 
an extremum of S3 (r ). A similar approximation can be used to estimate the 
radius of a thick- walled bubble (i.e., outside the range of validity of the thin 
wall approximation) |I3| . However, as was pointed out in Ref. [§], one must 
be very careful when evaluating S3 in this way. Since the solution of Eq. fllCf) 
is a maximum of the action, the corresponding value of S3 will be higher than 
the action of any approximate solution. As a consequence one can strongly 
overestimate the tunneling probability, leading to a sooner completion of the 
phase transition. 

For the potential (0), Eq. ( |T0"D can be solved numerically and then inte- 
grated to obtain the action S3, which can be expressed as a function of a 
dimensionless parameter a = ( T z j , 

| = 13.72^ a 3 / 2 /(«) . (11) 

Dine et al. |§ have found a useful analytical approximation to their numerical 
result for the function / (a), 



a / 2.4 0.26 

- H 1 

4 V 1-a {I -a) 



^Hi + tU + t^ + t^ti) , (12) 



with an accuracy of about 2% in the interval^ < a < 1. Using the same 
procedure we have calculated the radius r of the nucleated bubble, the 
thickness l w of its wall, and the value of the Higgs field inside it, <fio, as 
functions of a. The calculation is described in the Appendix. Notice that 
the dimensionless function p (T) defined by Eqs. and (^) has also a simple 
expression as a function of a, 



p{a) = - l + Wl-_ a . (13) 




Figure |l] illustrates the behavior of r (e) and l w (e). Near the critical 
temperature (e = 0) the bubbles that nucleate are thin- walled (l w r ). 
This is because the energy difference between the two minima of the potential 
is much smaller than the energy barrier between them. In such the 
radius must be very large, so that the negative potential energy inside the 
volume of the bubble is able to compete with the positive surface energy. 
As e increases, the wall width becomes of the same order of the radius and 



1 It is easy to see from Eq. (||) that a (Tc) = 1, a (To) = 0, and a ~ 1 — e in the whole 
interval. 
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the thin wall approximation breaks down. For e — > 1 the size of the critical 
bubbles diverges. This is because the width of the bubble wall is associated 
to the correlation length, l w ~ £ (T) ~ m (T) _1 , with m (T) = 2D (T 2 - T 2 ). 
This divergence is easily understood by looking at the shape of the bubbles. 
Figure |^ shows the bubble profile at a temperature near T c (e = 0.1), a 
temperature near T (e = 0.9), and at an intermediate temperature (e = 0.5). 
Note that the value of the Higgs field inside the bubble decreases with e. In 




1 ' ' ' 1 ' ' ' 1 ' ' ' 1 1 ' ' 1 ' ' 1 1 

0.0 0.2 0.4 0.6 0.8 1.0 

£ 

Figure 1: Bubble radius and wall thickness of the nucleated bubbles vs the 
relative messure of the temperature, e, for D = 1, E = 0.06, and A = 2E. 




r(T ') 

Figure 2: Shape of the critical bubble at different temperatures. 

Fig. |H we have plotted this value together with the minimum of the potential 
v (T) for comparison. At e = the value of 0o coincides with that of v (T c ). 
However, as the temperature decreases from the critical one, this value moves 
away from the minimum of V (4>) and goes to zero. This behavior is due 
the decreasing of the height of the barrier and the increasing in the energy 
difference between minima. In this case the value of (x) inside the bubble 
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needs not anymore to be exactly at the minimum of the potential to get 
enough volume energy. Moreover, as the barrier between minima disappears, 
it becomes easier to form a large bubble with a small value of inside it. 
Therefore, near the temperature T , bubbles with r — > oo and O — > will 
be formed. 




Figure 3: The minimum of the potential v (T) and the value of inside the 
nucleated bubble, 0o (T) = (r = 0). 

Consequently, if for any reason the transition does not occur until the 
temperature is close to T , then it will tend to proceed by a homogeneous 
growing of the order parameter 0, instead of by the nucleation and expan- 
sion of bubbles. This is not surprising since, without a barrier, one expects 
the field to "roll down" homogeneously towards the minimum of the effective 
potential ||18|| . There would still be some baryon number generation in such a 
scenario, due to the coherent variation of the CP violating field (t), which 
would produce an effective chemical potential for baryon number. However, 
this would be a kind of local baryogenesis |IJ mechanism, and hence it would 
lack the enhancement due to diffusion, which characterizes the usual non- 
local electroweak baryogenesis. Furthermore, after rolling down the field will 
start to oscillate around the minimum of the potential until all its energy 
is transferred to the plasma. This oscillations can further spoil the gener- 
ated baryon asymmetry, depending on their amplitude. This effect would 
be avoided in a phase transition at a higher scale, where the motion of 
would be strongly damped by the expansion of the Universe apart from the 
damping due to the plasma. Nevertheless, we shall see in section ^ that the 
phase transition is unlikely to occur beyond e = 0.5, so the value of inside 
the bubbles at the moment of nucleation will not be far from v (T) . 

We also see from Fig. |I| that the wall thickness is almost constant up to 
that value of e, so we can anticipate that the width l w , (which enters the 
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generated baryon number), will depend weakly on the temperature of the 
phase transition. 



3 The evolution of bubbles 

After a bubble is formed, it starts to grow due to the pressure difference 
between the phases separated by the bubble wall, AV — V (v, T) — V (0, T). 
After a short acceleration stage, the wall reaches a terminal velocity which 
depends on its interaction with the plasma. In this section we will analyze 
the evolution of bubbles from the moment they are nucleated until they fill 
the Universe 



3.1 Bubble wall propagation 

If the time it takes a bubble wall to accelerate up to terminal velocity is much 
smaller than the time it propagates before colliding with other bubbles, then 
the initial acceleration stage can be ignored. In order to check that this is the 
case, we can estimate the time needed for the wall to reach terminal velocity 
from the equation of motion for <fr in the plasma, Eq. (|i~7D below. If the field 
is strongly damped the wall will reach terminal velocity faster than without 
damping, so using Eq. flT7|) without the damping term one can estimate the 
maximum acceleration time ||. In such the wall reaches the speed of 

light in a time r ~ v (T) 2 / (l w V (v (T) , T)). For the values of the parameters 
of V (4>, T) that we will consider, this gives r ~ (10 2 — 10 4 ) T~ l . 

In order to compare this time with the typical times associated with 
the phase transition, we need the relation between the temperature of the 
Universe and time. At the radiation-dominated era it is given by 

*=^. d4) 

where Mpi = 1.22 x 10 19 Ge^ is the Plank mass, and £ is a factor that depends 
on the effective degrees of freedom of the plasma. Near the temperature of 
the electroweak phase transition, £ ~ 1/34. The time elapsed between T c 
and T , for instance, is 

At Tc -T ^ ^ J-c ~ 10 1 , (15) 

which is many orders of magnitude larger than r. However, to check if the 
short acceleration assumption is correct we must compare r with the time in 
which a typical bubble propagates until all space is filled by broken symmetry 
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bubbles. This time is of the order of the duration of the phase transition Att, 
which is less than At Tc _ To . We will see in section [| that AT t ~ 1CT 3 (T c — T ), 
so At t ~ 10 10 T _1 . The assumption is thus valid. 

There is another transitory stage, associated with the growing of in the 
interior of the bubble. As seen in Fig. |3], bubbles nucleate with a value of the 
field inside them, 0o (T), that is less than that of the minimum of the poten- 
tial v (T). This fact can be important for the bubble wall propagation, since 
the pressure difference V ((f>o, T) — V (0, T) that pushes the wall would then 
be less than V (v, T) — V (0, T). Again, we can use Eq. (|17D , this time dis- 
carding all spatial derivatives, to calculate the time in which (r = 0) grows 
from (J)q to v . Neglecting the damping term we obtain r ~ ^A(f)/V'(<p). For 
A0 ~ 0o — v/2 this gives r ~ 10T -1 . With a damping term it will take the 
field a longer time to get to v. Assuming a typical damping rj ~ 100 (see 
below), we obtain r ~ A^rj/V ((f)) ~ (10 4 — 10 5 ) T , which is again much 
less than At t . 

We can therefore safely assume that (r = 0) — v (T), and that bubble 
walls move with constant velocity through the plasma. The friction on the 
bubble wall depends on the departure from equilibrium of the particle distri- 
butions near the wall || [19|, p0| . Particles which interact with the Higgs field 
in such a way that their masses change across the bubble wall profile will feel 
a force as they are caught by the propagating wall. Integration of these forces 
weighted with the particle distributions gives the total force exerted on the 
wall by the plasma. When equilibrium distributions appropriate for the local 
value of are considered, this force only accounts for the finite temperature 
part of AV, and does not depend on the wall velocity. Nevertheless, a depar- 
ture from equilibrium is in fact produced by the motion of the wall. Particles 
that, for energetic reasons, cannot penetrate the wall, are reflected in front 
of it. These particles slow down the wall by giving momentum to it, and this 
effect becomes more important the faster the wall moves. Expanding the 
non-equilibrium particle distributions to first order in the wall velocity, the 
force gets a contribution proportional to the wall velocity ||. The terminal 
velocity is thus obtained by equating this frictional force with AV. 

The propagation of the bubble wall can also be affected by hydrodynamics 
(see for instance || [22]). Although hydrodynamic considerations are 
important for large velocities, if the wall velocity is small the only important 
effect is a global reheating due to the latent heat released by the expanding 
bubbles 0. Indeed, when the velocity of the bubble wall is less than the 
speed of sound in the relativistic plasma c s = a/1/3, it propagates as a 
deflagration front. This means that a shock front precedes the wall, with a 
velocity v s h > c s . Hence, the latent heat is transmitted away from the wall. 
If the wall velocity is small enough, the latent heat has time to distribute 



11 



throughout space, causing a homogeneous reheating of the plasma as the wall 
propagates. 

There exist other effects that can influence the growth of bubbles. An 
example is that of strong magnetic fields that can be formed before [E3j or 



during |24|] the electroweak phase transition. The presence of these fields 
changes the free energy difference between phases, thus affecting the temper- 
ature and strength of the first order phase transition ||25|| . In the case that the 
magnetic fields are generated during the electroweak phase transition, they 
can delay the completion of the latter by strongly affecting the evolution of 
bubbles 



3.2 The wall velocity 

The equation of motion for the bubble wall in the hot plasma can be derived 



by energy conservation considerations H, 19, GO 



d<\> J (2ttY2E 1 

where Sfi is the deviation from the equilibrium distribution function for 
particle species i in the heat bath. The last term can be expressed || |21 



as 



a frictional damping term proportional to dcj)/dt. Due to Lorentz invariance 
this term must be of the form u^d^cf), where is the four-velocity of the 
plasma. Eq. flTop then becomes 

8^6 + V (6) + (rjT) u^S = , (17) 

where rj is a dimensionless damping coefficient that depends on the viscosity 
of the medium, and is obtained by calculating the deviations Sfi near the 
wall. 



The calculation of the friction is usually treated with kinetic theory M, [19 



20] , giving for the Standard Model roughly r\ ~ 1. However, it was noticed 
]2"6] | that infrared boson excitations which are treated classically [|27|| , increase 
the value of r\ by an order of magnitude. In the MSSM there is an additional 
enhancement of the friction due to the larger particle content of the model. 



In Ref. [^] the contribution of a light right-handed stop has been calculated, 
which added to the contributions of tops and W bosons give a lower bound 
on the friction coefficient of rj ~ 70. 

We will not take into account other effects such as the influence of mag- 
netic fields on the bubble growth. Instead, we will consider values of the 
friction beyond rj ~ 100 to allow for the possibility of additional damping, 
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and we will assume that that kind of effects can be parametrized in this way. 
Accordingly, we will let rj vary in the range 1 < rj < 1000. 

The nucleated bubbles grow rapidly to a size large enough so that we can 
go to 1+1 dimensions. Boosting to a frame that moves with the wall and 
assuming stationary motion in the x direction we then have 

^V»)-^, (18) 
dx A ax 



where 7 = 1/y/l — v 2 . Multiplying both sides by d(f)/dx and integrating over 
— 00 < x < 00, we obtain a formula for the wall velocity, 

1 AV 

Vwl = -= , (19) 

r/T a 

where AV (T) is the free energy difference between the two phases and a (T) 
is the surface tension of the wall, 

+ °° / rlM\ 2 



° {T)= L\f x ) dx - (20) 

We have assumed here that temperature is constant across the wall. This 
is right if the wall velocity is small enough that the latent heat it releases 
has time to be uniformly distributed throughout space || . If we use the well 
known kink approximation for the wall profile (see for instance |3(J ) 



1 + tanh 



(21) 



then the surface energy ( f2~0|) is a = v 2 /3l w , and the wall velocity is given by 



ivy = - 2 • 22 

We see from Eq. ( ]2"2] ) that v w increases with l w . This is because the 
thicker the bubble wall, the slower will be the variation of <p as the wall 
sweeps through a given point in the plasma, and consequently, the smaller 
the disturbance of the plasma near the wall. Hence, thicker walls experience a 
smaller friction force from the medium. Since l w diverges at T , one can infer 
that the dependence of the velocity with the wall width gives an important 
variation of v w with temperature. In the literature, however, l w is usually 
evaluated near the critical temperature. For instance, it is often roughly 
approximated by the correlation length m (T c ) _1 [[3], |26|, or it is taken from 
the wall profile at the temperature at which S 3 (T) /T ~ 130-140 0, §, |, [28 . 
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This gives l w ~ (10 — 40) T -1 for the electroweak bubble. We can see from 
Fig. |I] that these approximations are good as long as the phase transition 
completes before the temperature gets close to T , since l w is almost constant 
until £ ~ 0.8. As we will see in the next section, the phase transition hardly 
occurs beyond that value of e. 

Figure [| shows the terminal velocity of the bubble wall as a function of 
e. In the interval < e < 1 it varies from v w — (since AV = at the 
critical temperature) to v w = 1 (since /„, — > oo for e — > 1). Note, however, 
that at the beginning of the interval the wall velocity grows rapidly from 
to some value proportional to r? -1 , then it becomes rather insensitive to e, 
and finally it grows again up to v w = 1. As we shall see, the electroweak 
phase transition generally occurs at some intermediate e so, for fixed rj, the 
wall velocity will not be very sensitive to the temperature of the transition. 
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Figure 4: Bubble wall velocity as a function of e, for three different values of 
the damping coefficient rj. 

3.3 Fraction of space occupied by bubbles 

Once a bubble has nucleated, its wall rapidly reaches the terminal velocity 
v w , which depends on temperature. Due to the large time-temperature re- 
lationship flI3|), v w varies slowly with time, so the evolution of the bubble 
radius can be approximated by 

r(t',t)=r (t')+v w (t')(t-t') (23) 

where r (£') is the radius of the bubble at the moment t' of its formation^. 
The number of bubbles created per unit volume between t' and t' + dt' 

2 By the epoch of the electroweak phase transition, the variation of length scales due to 
the expansion of the Universe can be neglected during the short period between T c and 
T . 
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is drib (f) = T (t') dt'. Then at time t the total volume of bubbles that 
have nucleated at previous times t c < t' < t is (per unit volume of space) 
f T (t 1 ) (47r/3) r (t 1 , t) 3 dt'. This formula, however, contains multiple count- 
ing of volume where bubbles overlap. A more exact expression was obtained 
in Ref. [29] for the fraction of space / (t) that is still in the false vacuum, 

/ (t) = exp {— Cv (f ) [r (f ) + v w (f ) (t - t')f • (24) 



The phase transition then completes when / (t) falls to zero. All the quanti- 
ties appearing in the integral are in fact functions of temperature, so the use 
of the time-temperature relation ([14]) will be needed in the computation. 

Some reheating will occur at the end of the transition due to collisions 
of bubble walls. However, for small wall velocities this effect should not 
be important. As we have already mentioned, the main consequence of the 
release of energy by slowly-expanding bubbles is a global reheating during the 
bubble expansion stage, which would manifest itself by changing the time- 
temperature relation. Taking into account this effect would require solving 
the two coupled equations for T (t) and / (t) , which is outside the scope of 
the present paper. Before going on, however, we would like to discuss the 
relevance of such an effect. As pointed out in Ref. 0, one can get an idea of 
how important the reheating is by comparing the latent heat 



T 



dV(<f>,T) 



DE 



2 



v,T c A 



n (25) 



dT 

with the energy needed to bring the plasma back up to T c from T t , 

p (T c ) - p (T t ) = a (T c 4 - , (26) 

where the factor a depends on the number of effectively massless degrees 
of freedom. At the time of the electroweak phase transition a ~ 35. Since 
T c — T t < T c — T , the energy density difference will be Ap < T c 4 . For 
SM values of D and E, the latent heat is smaller than Ap by one or two 
orders of magnitude, depending on the value of A. For MSSM values of the 
parameters, instead, L can be of the same order of magnitude of Ap. In such 
a case the latent heat released by the bubbles could reheat the plasma back 
up to temperatures near T c . So we must remark that this effect might be 
important for baryogenesis. Due to the small velocity of the bubble walls, 
the gradually released latent heat will prevent the temperature of the bath to 
decrease, rather than reheating the plasma at the end of the transition. If the 
temperature stays close to T c , the pressure difference Ay across the bubble 
walls will be small, and hence the wall velocity can decrease considerably 
from the current estimations 131]. 
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4 The temperature of the phase transition 



The progress of the phase transition is determined by the fraction of space 
/ (T) that is still in the symmetric phase, given by Eq. (|24[) . In the next 
section we will calculate the temperature T t of the transition and its duration 
AT t for different values of the friction. However, since the exact evolution 
of the bubble wall velocity is not known, it will be convenient to begin our 
analysis by putting reliable constraints on T t inside the interval T — T c . 

We will consider two sets of values for the parameters E and D of the 
potential: those corresponding to the SM, E ~ 6 x 1CT 3 , D ~ 0.2; and 
MSSM-like values, which for defmiteness we set E ~ 6 x 10~ 2 , D ~ 1. In 
both cases we will take X/E = 2 or 3, which correspond respectively to 
v (T c ) /T c = 1 or v (T ) /T = 1. Although the SM with these values of the 
parameters is not a realistic model, for the sake of comparison it is useful to 
consider such a kind of parameters besides those corresponding to the MSSM, 
in order to simulate a different phase transition with the same simple form 
of the effective potential. Regarding the MSSM-like values, it is worthy to 
remind that our model Eq. (|2J) corresponds to the limit of only one light 
Higgs boson of the MSSM, and the values of the parameters correspond to 
the light right-handed stop scenario. 



4.1 An upper bound 

The temperature T/v at which nucleation of bubbles begins is defined as 
that at which the first bubble is formed inside a causal volume In the 
radiation dominated era the horizon size scales like du = 2t, where the age 
of the Universe t is given by Eq. (|14]). Hence the size of a causal volume is 
Vh ~ 8£ 3 M Pi /T 6 , and the probability of nucleating a bubble inside it is 

ft rT irp 

P (T) = / YV H dt = (2£M PZ ) 4 / e~^ T — , (27) 

Jt c JT C J- 

where we have used Eqs. (|8]) and (|TJ) in the last step. The temperature T/v 
is thus determined by the condition P (T/v) ~ 1. This can be solved either 
by an analytical approximation [13] or numerically, and gives 



S 3 (Tv) /Tv ~ 135 . (28) 

Although T/v is often assumed to be the temperature of the transition, it 
only corresponds to the onset of nucleation, so it gives an upper bound for 



16 



Once bubbles begin to nucleate, however, it still takes them some time 
to fill all space, since bubble walls cannot move faster than light. A less 
conservative bound T max < T/v can thus be obtained by setting v w = 1 instead 
of using Eq.(^) for v w in Eq. ([24|). In this limit the expansion of bubbles is 
so fast that we can neglect their initial radius r . Then the condition that 
the fraction / (T) falls to zero can be expressed as 

f {2iM Pl f £° e**™ {T - Tf ^ ~ 1 , (29) 

giving for the upper bound T max 

S 3 (T max ) /T max ~ 105 - 109. (30) 

The small uncertainty is due to the variation of the parameters of the poten- 
tial as we take SM-like or MSSM-like values and the ratio X/E in the range 
2-3. 



4.2 A lower bound 

It is easy to see from Eqs. (|TT|) and (|12"D that the instanton action S3 vanishes 
at To. Hence, there is no exponential suppression to the nucleation rate T at 
that temperature (this is a consequence of the disappearance of the barrier 
of V (4>) at Tq). At temperatures near T the nucleation rate is thus F ~ T 4 , 
which is extremely large. Indeed, T = T 4 means that in a time interval 
t ~ T _1 a bubble forms in every volume V ~ T~ 3 . This time r is many 
orders of magnitude less than the time elapsed between T c and T , Eq. (|T5|). 
In addition, the size of the nucleated bubbles is larger than the volume per 
bubble V, since the radius of the electroweak bubble is always larger than a 
few T _1 . So all space must be filled with bubbles before the Universe cools 
down to some temperature T min > T , when the exponential factor in Eq. (|3|) 
is still much less than 1. 

We can calculate the lower bound T min by setting v w = in Eq. (|24]), 

/ r(T)r 3 (T) — ~1. (31) 

" ^min 

This corresponds to the (unrealistic) limit of a nucleation-dominated first- 
order phase transition. The opposite limit is that of a bubble expansion 
dominated phase transition, in which the nucleated bubbles grow so rapidly 
that the original size is negligible, and the transition is completed when the 



nucleation rate is still small fT3]. This is the case of the upper bound T max we 
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have just estimated. As we shall see, the electroweak phase transition is closer 
to an expansion dominated transition than to a nucleation dominated one. 
The temperature T min determined by Eq. fl3Ip is given by S3 (T min ) /T min ~ 
35 - 38. 

In terms of the relative variable e, the bounds on the temperature of the 
transition are expressed as 0.31 < e t < 0.52 for SM values of E and D, and 
0.15 < Et < 0.31 for MSSM values of the parameters (setting for defmiteness 
A = 2E in both cases). 

5 Numerical results 

In the last section we have determined the temperatures T max and T min be- 
tween which the phase transition takes place. Depending on the value of the 
friction parameter 77, the temperature of the phase transition will be closer 
to one of these bounds. 

5.1 The development of the phase transition 

The fraction of space that is still in the symmetric phase is obtained by in- 
serting Eq. ( |2"2"D into Eq. (^) and solving numerically for f(T). We have 
plotted / (e) in Fig. || for a given choice of parameters, in order to illustrate 
the evolution of the phase transition. Note that, once the fraction of volume 
occupied by bubbles becomes appreciable, the phase transition completes in 
a tiny temperature interval. That is why we speak of the "temperature of the 
transition" rather than of a "temperature interval" . No relevant quantity will 
change appreciably in this interval. This is due to the large time-temperature 
relation, which governs the nucleation and expansion of bubbles. The du- 
ration of the phase transition varies very little with the friction. It is given 
by 

AT 

* -A £f ~5x 10~ 3 (32) 



T c -T 

for values of 77 as large as 10 5 as well as for no friction at all. 

With no friction, the phase transition occurs at T max . When friction is 
turned on, the transition is retarded, as shown in Fig. || Note that the 
dependence of e t with rj is almost logarithmic, so e t increases more slowly for 
large 77. As a consequence, T t will be closer to T max rather than to T min for 
not extremely high values of 77. As discussed in section ||, 77 should not be 
out of the range 1 — 1000. 

In Fig. ^we have plotted the value of S3/T at the moment of the transi- 
tion, as a function of the damping coefficient. The solid line corresponds to 
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Figure 5: Fraction of space in false vacuum as a function of e, for rj = 100, 
D = 0.2, E = 0.006 and A = IE. 
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Figure 6: The temperature of the transition as a function of the friction rj. 
The upper curve corresponds D = 0.2, E = 6 x 10~ 3 , and A = IE. The lower 
one is for D — 1, E — 6 x 10 -2 . The temperatures corresponding to upper 
and lower bounds are indicated with dashed lines. 

MSSM values of the parameters. For likely values of the friction, rj ~ 100, 
the temperature of the transition is given by 53 (T t ) /T t ~ 90. The dashed 
line indicates the variation in S3/T when considering different parameters of 
the effective potential, as explained in the previous section. Since the ambi- 
guity in 5*3 (T t ) /T t produced by such a variation is of only a 5%, we believe 
that Fig. not only describes qualitatively the general case, but it would 
also give a good approximation for the determination of the temperature of 
the electroweak phase transition within different models. 
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Figure 7: The value of the three dimensional action S3 at the temperature of 
the transition, for D = 1, E = 0.06 and A = 2E. The dashed line shows the 
maximum shift of the curve as D, E and A are varied as indicated in Sec. |1. 



5.2 Implications for baryogenesis 

The amount of baryon number generated at the electroweak phase transition 
depends on the shape of bubbles. We have seen that the relevant parameters, 
namely, the width and velocity of the bubble wall and the value of the Higgs 
field in the broken symmetry phase, vary significantly between the tempera- 
tures T c and To, and hence the importance of determining their values at the 
intermediate temperature T t at which the phase transition occurs. 

Regarding the wall velocity, we have shown in section |3| that, although 
v w takes all the possible values < v w < 1 in the interval < e < 1, it 
changes rapidly at the beginning and at the end of that interval, but remains 
of the same order of magnitude for intermediate values of e, roughly between 
e = 0.1 and e = 0.8 (see Fig. |]). Since the bounds on T t obtained in section | 
fall in this region of slow variation of v w , we see that the wall velocity depends 
weakly on the actual temperature of the transition, although it is sensitive 
to the friction, as shown in Fig. |j. Therefore, rough approximations of the 
temperature of the transition, such as the usual T t — Tn, do not introduce 
significant errors in the estimation of v w . The same can be said of the wall 
thickness, as inferred from Fig. [l]. 

Typical values of these quantities are found by taking for instance r\ = 70 
and MSSM-like parameters for AV^. Then we obtain l w (T t ) ~ 16T -1 , which 
inserted into Eq. (|22|) gives v w (T t ) ~ 3 x 10~ 3 . This is in agreement with the 



value obtained in Ref. [2q|, and supports an MSSM scenario for electroweak 
baryogenesis. 

With respect to the Higgs mean value v (T), although it varies by a 50% 
between T c and T , the variation between the critical temperature and T t 
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will be smaller. As shown in Fig || the temperature of the transition stays 
close to the upper bound (lower bound on e t ), even for large values of r\. 
We can use Eq. (|TB|) to estimate by which amount the order parameter can 
change. Setting for defmiteness A = 2E, like in Fig. || we see that for rj ~ 100 
the upper curve in the figure gives e t ~ 0.35. Since v (T) /T = (E/X)p(a) 
and a ~ 1 — e, we readily obtain v (T t ) /T t > 1.2. Similarly, for the lower 
curve we obtain v (T t ) /T t > 1.1. Comparing with v (T c ) /T c = 1, we find a 
difference of at least a 10 — 20% between the critical temperature and the 
actual temperature of the transition. This means that the use of Eq. ([!]) 
can result in an upper bound on the Higgs mass that is by that percentage 
more severe than the correct constraint. A more precise bound for A can be 
obtained by using Eq. (0) recursively. For instance, after computing Et with 
the assumption A = 2E and calculating the constraint A = p (Et) E, this new 
value of A can be used to recalculate e t and obtain a higher order correction 
to the bound. This may be useful when considering a more concrete model 
for the phase transition. 

The error in the determination of v (T t ) /T t diminishes considerably when 
using the temperature Tjv of the onset of nucleation instead of T c . For the 
values of the parameters that we have considered, v (T) /T increases by less 
than a 5% between T N and T t . Hence, regarding the constraints that result by 
avoiding the washout of the baryon asymmetry, we again find that T t ~ Tjy 
is a good approximation in the calculations of electroweak baryogenesis. 

Before concluding, we would like to discuss the behavior of another bubble 
parameter that is relevant for baryogenesis in the MSSM, that is, the relative 
variation of the two Higgs fields across the wall, A/3. The CP-violating force 
on the quarks depends on this variation, and thus the final baryon asym- 
metry results proportional to A/3. In order to investigate the temperature 
dependence of this parameter, however, we would need to consider the full 
two-Higgs model. The computation of the bubble profile with more than 
one scalar field is a difficult task (see for instance ||| |30|). Nevertheless, our 
results for the temperature of the transition can be used to verify the ap- 
proximations employed in existing calculations of Af3. Moreno et.al. B, for 
example, have calculated the Higgs profile of the bubble in the light stop sce- 
nario. They performed the calculation at the temperature we have called Tjy, 
and also at T N + 0AGeV and T N — OAGeV, and the variation they obtained of 
Aft along the bubble wall is very similar among these temperatures. Trans- 
lating that temperature interval to the variable e by Ae ~ AT/ (T c — To), 
where AT = 0.8GeV and T c — T ~ 3GeV, gives Ae ~ 0.2. So the conclusion 
is that the parameter Af3 does not change significantly in an ample interval 
around En, which, according to our results, includes the possible values of £ t . 
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6 Conclusions and discussion 



It is well know that phase transition dynamics plays an important role in elec- 
troweak baryogenesis. Although the production of the baryonic asymmetry 
of the Universe at the electroweak phase transition is an exciting possibility, 
it is clear that a complete picture of the phase transition is needed to check 
quantitatively this possibility. This includes the detailed investigation of the 
evolution of bubbles, as well as the accurate calculation of the temperature 
of the transition. In this paper we have investigated both of this problems, 
which are certainly related with each other. 

We have determined the upper and lower temperature limits T max and 
T min inside the interval T — T c , between which the phase transition must 
take place regardless of the velocity of expansion of bubbles. Our main 
conclusion in this respect is that, for the electroweak phase transition, the 
lower bound on T t avoids the possibility that the phase transition occurs at a 
temperature near T . This is good for baryogenesis, since in such a case the 
transition would proceed by an almost homogeneous growing of the Higgs 
field, what would lead to the generation of a smaller baryon abundance. 

When the evolution of the bubble wall is included, the temperature of 
the phase transition varies between the limits T max and T min as the friction 
parameter 77 goes from to 00. However, it turns out that, due to the 
slow decreasing of temperature with time, the temperature T t stays close 
to T max for not extremely high values of rj. According to our numerical 
computations, for likely values of the viscosity of the plasma the temperature 
of the electroweak phase transition is given by S3 (T t ) /T t < 90. 

On the other hand, the values of the wall thickness l w (T) and wall velocity 
v w (T) are rather insensitive to temperature in the range T m i n — T max , so they 
depend weakly on the actual temperature of the phase transition. Regarding 
the Higgs mean value v (T), although baryogenesis is very sensitive to this 
parameter, and hence the uncertainty between T min and T max would still 
be considerable, the result that the phase transition takes place near T max 
eliminates this uncertainty. Furthermore, the temperature T max is in general 
close to that of the onset of nucleation. Hence we conclude that l w (Tt), 
v w (T t ), and v(T t ), which are relevant parameters for baryogenesis, do not 
change significantly from those obtained with the usual estimation of the 
temperature of the transition, i.e., with the condition S3 (T) /T ~ 135. We 
have also shown that the same applies to other quantities, which we have 
not computed, such as the parameter Af3, which is relevant for electroweak 
baryogenesis in the context of the MSSM. 

Regarding further development of this analysis, a detailed investigation 
of the effect of reheating due to the latent heat released by the expanding 
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bubbles would be necessary |JT|. As we have discussed in section f|, such an 
investigation would be especially important in the case of the MSSM, where 
the main consequence of this effect would be to prevent the wall velocity 
to reach the value ~ 10~ 3 — 10~ 2 that is needed for a satisfactory result of 
electroweak baryogenesis. 
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A The shape of bubbles 

In this appendix we give some interesting details of the calculation of the 
bubble radius, the wall width, and the value of the Higgs field inside the 
bubble at the moment of formation. 

The shape of the bubble is given by Eq. ([U]) with potential (§), 

^ + 1 d ± = 2D (T 2 - T 2 ) - 3£T0 2 + A0 3 . (33) 
dr 2 r dr v 

Follo wing Dine et al. f§, we define = [D (T 2 - T 2 ) / (ET)] $ and r = 
R/^2D (T 2 -T 2 ), so that Eq. (H) becomes 



rf 2 $ 2# ^ 3^ 2 1 ^ 3 

dR> + RlR = ^t + 2 a$ ' (34) 

with a = \D (T 2 — T 2 ) / (ET) 2 . We have solved numerically this equation 
with the usual overshooting-undershooting method, to obtain the configu- 
ration $ (R) for different values of a. The graphics of Fig. |2| for 0/T as a 
function of the dimensionless radius Tr are then given by 




(35) 



The expression ( |TT| ) for S 3 (T) /T is then obtained by an appropriate inte- 
gration of the solution over d?X = 4nR 2 dR. 

We define the radius Rq (a) of the configuration as that at which $ (R) 
falls to $o/2, where $ ( a ) = $ (-R = 0). Similarly, we calculate the radii Ri 
and R 2 corresponding to 0.9$ an d O.l'I'o respectively, and define the wall 
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width as L w — R 2 — R\. This quantities are related to the physical radius r 
and width l w of Fig. [I] by 

Tr = y^J-R ( a ) , Tl w = ^Sl^-L w (a) . (36) 
L \/a hj J a. 



The value of the field inside the bubble, plotted in Fig. [| is defined as 
0o = <t> (r = 0) and is related to the dimensionless $ by 

^ = f «*.(«), (37) 
to be compared with the minimum v (T) of the potential, given by 

^ = !*>(«)■ (*» 

The divergence of r and i to at a = is apparent in Eq. (|36|) . The 
divergence of r at a = 1, instead, is hidden in the function R (a). The 
wall width is finite at a = 1, but due to the divergence of the radius the 
shape of the bubble wall cannot be calculated using the numerical method. 
Interestingly, in this limit L w can be obtained analytically, what is useful 
to check that the numerical curve goes to the right limit for a — > 1: Since 
R — > oo the second term disappears from Eq. fl3H|) . Then multiplying by ^ 
and integrating over dR gives for a = 1 

Integrating again across the wall we obtain 

'$ -2 



L M = In 



0.1*o 

(40) 

0.9*o 



Similarly we can see that $ (R) goes from 2 to across the wall, so in this 
case $o = 2, and Eq. ( |4"0"D gives L w = 2 In 9. We can then verify that in 
Fig. [I], with A = 2E and E = 0.06, the wall width has the correct value at 
£ = 0. 
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